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Abstract 

We present an exact formula for the computation of the interference rate of radi- 
ation in the case of two charges revolving with constant angular velocity at opposite 
ends of a diameter in a fixed circle. The formula is valid for arbitrary velocities 
of the charges, and can be easily studied by numerical methods, even for veloci- 
ties very close to the velocity of light. For ultrarelativistic motion, the interference 
rate of radiation behaves as ln(l — w-^/c^)"^/-^, which contrasts with the behavior 
(1 — for the rate of radiation for one charge in circular motion. This is 

the first exact calculation for the interference rate of radiation of two relativistic 
charges, and it is useful in connection with the old controversy about the correctness 
of the Lorentz-Dirac equations of motion for more than one charge. 
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I Introduction 



In the case of a one charge in arbitrary motion, the total rate of radiation emitted 
at time t is given by the Larmor's formula, which reads 

^ = ^^/{(/3)^ + f(/3-/3)^}, (1-1) 

where v is the charge velocity, c is the velocity of light, /3 = v/c, (3 = dp/dt, 
and 7 = (1 — The simplicity of this formula, where the variables (3 and /3 

are evaluated at the same time as the total rate of radiation, reflects some special 
properties of the point charge field. There are basically two derivations of the 
Larmor's formula ( |1.1| ). In one of them ^, the total rate of radiation is first 
calculated in the Lorentz rest frame of the charge, and equation (^]^) is obtained 
using covariance arguments under Lorentz transformation. The other derivation is 
carried out directly in Minkowski space, using the fact that the radiation of a point 
charge can be characterized locally [§, |5|. In particular, because of this property, 
it is not necessary to go very far from the charge in order to calculate the total rate 
or radiation. 

In constrast with the one charge formula for the total rate of radiation of 

two charges under arbitrary motion is not known yet. Furthermore, does not exist in 
the literature an exact formula for the total rate of radiation for any special type of 
motion of two charges. The source of the technical difficulty comes, of course, from 
the fact that the fields of the charges are explicit functions of the retarded times 
instead of the laboratory time. This point can also be illustrated by referring to the 
already mentioned derivations of ( [1.11 ). Thus, as emphazised by Landau & Lifschitz 

for two charges there is generally no system of reference in which both charges 
are at rest simultaneously. And, on the other hand, the local characterization of 
radiation does not work for more than one charge Q. 

Due to the superposition principle and the quadratic nature of the Poynting's 
vector, the energy flux for two charges contains three terms. Two of them are asso- 
ciated with the field of each individual charge separately. The third one corresponds 
to an interference term that mixes the fields of both charges. In the evaluation of 
the total rate of radiation, the first two terms give rise to a Larmor formula for each 
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charge; so the real problem is the calculation of the energy flux, across the surface 
of a sphere of very large radius, of the interference term. 

There have been only a few attemps to calculate the total rate of radiation as- 
sociated with the interference term. Huschilt and Baylis |^ studied this radiation 
in the case of two identical charged particles, that are moving in a straight line in 
head-on collision. This calculation is carried out with the help of the equation of 
motion for the charges, and in addition it involves some kind of non-relativistic ap- 
proximations. Aguirregabiria and Bel studied the radiation of two charges using 
a covariant formulation. These authors calculated an integral of the interference 
field over a circle, for a rather general motion of the charges. From this result, 
they elaborated a formalism, with the help of the equations of motion and some 
additional assumptions, in order to obtain succesive approximations for the total 
rate of radiation. The interference radiation has been studied also by Hojman et al 
[P, by means of a covariant formalism. However, these authors also introduce some 
kind of non-relativistic approximations. 

In this paper we present an exact formula for the total rate of radiation in the 
case of two charges moving in a plane at the opposite ends of a diameter, revolving 
at constant angular velocity in a fixed circular orbit. The main motivation for this 
calculation is that it helps to solve an old controversy about the correctness of the 
Lorentz-Dirac equations of motion for more than one charge [|lO|. In fact, we have 
recently showed that, with appropriate external fields, the Lorentz-Dirac equations 
describes the circular motion under consideration in the case of two particles of equal 



charge and mass [11]. Then, if we know the total rate of radiation for this motion of 
the charges, we can check the consistency of the Lorentz-Dirac equations with the 
energy conservation law. We carried out such an analysis in |Jl|, and showed that 
this type of circular motion allows us to see in a manifiest way the inconsistency of 
the Lorentz-Dirac equations for more than one charge with the energy conservation 
law. 

The calculation is carried out directly in the laboratory frame. This is more 
appropiate than the covariant techniques, since the total rate of radiation can be 
clearly visualized from a physical point of view, and the result is free of any ambigui- 
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ties whatsoever. Our formula is an integral expression for the total rate of radiation. 
Unfortunately, the integral is too complicated for closed analytical evaluation. Nev- 
ertheless, it can be easily studied using numerical methods, even for velocities near 
the velocity of light. Furthermore, the integral can be evaluated approximately for 
ultrarelativistic motion of the charges. We find that the rate of interference radia- 
tion grows as (/9^/4) In 771 when /? tends to one. This result differs strongly with 
the behavior of the Larmor term of each charge, which behaves as /3^7^ when /3 is 
near one. Thus, for example, if we consider electrons of 500 Mev, the interference 
radiation is completely negligible in comparison with the Larmor term, since this 
latter is of the order of 10^^, while the interference term is near 1. For low velocities, 
however, it is known that the Larmor and the interference terms are comparable 

In section II we evaluate the energy flux across the surface of a sphere centered 
at the orbit center of the two charges, and show that this flux is independent of the 
time at which it is evaluated, as well as of the radius of the sphere. In section III 
we present a power series expansion in /3 of the total rate of radiation up to /3^. 
In section IV we derive an exact formula for the interference radiation term. In 
section V we present an analytical approximation of the exact formula for the case 
of ultra-relativistic motion. 

II The interference rate of radiation 

In the following we will be concerned only with the radiation of two charges moving 
in a plane at opposite ends of a diameter, revolving at constant angular velocity to, 
in a fixed circular orbit of radius a. In this case, because of the symmetries of the 
motion, it is possible to identify without any ambiguities the total rate of radiation. 
With this purpose in mind we work directly in the laboratory frame, since this 
allows us to have a clear physical picture of the radiation. 

FIGURE 1 

As it is shown in figure 1, our coordinate system is such that its origin coincides 
with the center of the orbit, and the X — Y plane is precisely the orbit plane. In this 
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figure we have drawn the positions of the charges at an arbitrary time t, and two 
spherical surfaces Si, and E2 centered at the origin of radii ri and r2 respectively, 
with r2 > ri > a, where a is the orbit radius. For a given time t, the electric and 
magnetic fields E and B change in a very complicated way from one point to another 
over the surface Si; this is because the retarded times of the two charges, which 
are in general different, change in a complicated way with the position over Si. In 
particular then, the Poynting vector S = (c/47r)E x B is a complicated function over 
Si, and so it is the total flux across Si at time t that we are interested in. Since 
the charges arc moving jointly at constant angular velocity, the position of them at 
different times looks the same with respect to the whole surface Si. This means 
that the energy flux across Si, cannot depend on time. In Section IV we present a 
rigourous proof of this property, using the explicit form of the electromagnetic field 
of the two charges. 

Now, if we denote by u{x, t) the energy density of the electromagnetic field, and 
by S the Poynting vector (c/47r)E x B; then, since the domain Q bounded by the two 
spherical surfaces Si and S2 of Fig. 1 is free of charges, the following conservation 
law holds in it: 

V ■S + du/dt = 0. (2.1) 

From this equation we obtain 

/ (S-f)dSi-/ {S-r)d^2 = -i f u{x,t)d^x, (2.2) 

where f = (sin ^ cos (y?, sin sin 99, cos 9) is the unit normal to Si and S2. But, 
because of the symmetries of the motion of the two charges, it is clear that the total 
energy contained in is independent of time. Therefore the integral 

^(S-i-)a!S, (2.3) 

over the surface of the sphere of radius r is not only independent of time, but it is 
also independent of the radius r. In particular then, the integral (2.3) represents 
the total rate of radiation that escapes to infinity; and it can be evaluated over 

the surface of any sphere of arbitrary radius r, with r > a. Thus, for this special 
type of motion of the two charges, a reminiscence of the local characterization of the 
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radiation still survives, in the sense that it is not necessary to go to infinity in order 
to calculate the total rate of radiation. 

We remark that the energy flux (2.3) is also independent of time and of the 
radius of S, in the case of only one charge in circular motion with constant velocity. 
Therefore, if we denote by Ei, Bi, and E2, B2 the electric and magnetic fields of 
the charges ei and 62 respectively, the integral 

/■TT f'2'n 

{cr'^/Ai:) / slnOde / (Ei x B2 + E2 x Bi) • f- ^99, (2.4) 
JO Jo 

is independent of time and of the radius r. The integral (2.4) represents physically 
the rate of radiation due to the interference of the fields of both charges, and it will 
be discussed in detail in section IV. 

Ill Power series representation for low veloc- 
ities 

The total rate of radiation of several charges can be in principle calculated by 
referring the retarded times of the charges, to the actual time t, by means of a 
power series expansion in 1/c. The electric intensity of a group of charges can be 
represented by the following formula when r goes to infinity |12|. 



rc 



d \ 1 d2 / 



(3.1) 

X f > X f , 



where f = i sin 9 cos (p + j sin 6 sin 99 + k cos 9, and denote the position and 
velocity of the charge at time t. The far magnetic fields is B = f x E. 
In our case we have 

ri(t) = iacoswt + jasina;t, foo\ 
^2{'t) = —lacoscat—jasmiot. 

For this motion we have computed the series (3.1) up to terms of power c~^^, and 
the result is given in appendix A. Using these formulae, the calculation of the 
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energy flux across the surface of a sphere of very large radius r is straightforward, 
but cumbersome. The result is the following 



dWrad 2 ef c 4 



dt 3 a 



{l + 2/?2 + 3/?^ + 4/?^ + 5/?^ + • • •} 



+^ fl + 2/32 + 3/3^ + 4/36 + 5/38 + .. .} 

3 J 



7 ' 945 ' 10395 



(3.3) 

The first two series correspond to power series expansions of the Larmor term of 
each charge. In fact, for a charge in circular motion with constant velocity, the 
Larmor formula (1.1) is reduced to 

|^/3^(l-/32)-2. (3.4) 
3 a"^ 

The series 

,3.5, 

of (3.3) represents, of course, the rate of radiation associated with the interference 
between the fields of both charges. The first two terms of (3.5) are already known 
[12|. As we will show, the series (3.5) gives the total rate of interference radiation 



with an error less than 4% for /3 < 0.2. Clearly, the error increases with /?, and a 
power series in (3 for the interference radiation is hopeless for f3 near 1, since a very 
large number of terms would be needed in this region. 

IV An exact formula 

We will consider now the energy flux of the interference term across the band be- 
tween the angles 6 and 6 + dO over the surface of the sphere of radius r at time t, 
that is, the contribution 



2-K 

(El X B2 + E2 X Bi) -Mt^, (4.1) 





Z7l 

/ 
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to the rate given in Eq. (2.4). The electric field Ei generated by the charge ei is 
given by the well-known Lienard-Wiechert formula [|l3| 
^(n,-/3i)(l-/32)' 



Ei(x,t) = ei 



c 

ret 



{(ni-/3i) x/3i} 



• (4.2) 

ret 



The corresponding magnetic induction Bi is 

Bi = nixEi. (4.3) 

In equation (4.2) and (4.3), ni is the unit vector that points from the retarded 
position ri(ti) of charge ei, associated with the point x and time t, to the detection 
point X. The vectors (3i and (3i have already been defined in connection with 
formula (1.1); but now they have to be evaluated at the retarded position ri(ti) of 
the charge ei. Moreover, Ri represents the distance between the detection point x 
and the retarded position ri(ti) of charge ei; and ki denotes the following positive 
number 

Ki = l-ni-/3i. (4.4) 

The electric field E2 and magnetic induction B2 of the charge 62 are given by (4.2) 
and (4.3) respectively, but where the quantities n2, /^2; ^2 and K2 are referred 
to the retarded time t2 of the charge 62 associated with the detection point x and 
time t. 

In what follows we will carry out our calculations in the coordinate system shown 
in Fig. 2; where the positions of charges ei and 62 are described by the vectors Ti(t) 
and r2(t) defined in Eqs. (3.2), and the detection point is x = irsin^cosy? + 
jrsin^sinc^ + krcos0. In Fig. 2 we have drawn the positions of the charges at 
three differents times, namely at t, ti and ^2; where t is the time at which we are 
going to calculate the fiux across the surface of the sphere of radius r > a; the time 
ti corresponds to the retarded time of charge ei associated with (x, t), and ^2 is the 
retarded time of charge 62 associated with (x,t). 



FIGURE 2 



Since the time interval t — ti that needs ei to go from its retarded position Bi 
to the actual position Ai is the same that takes the light for travel from Bi to P, 
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we have 

t-ti = {a/c)C^ {l + ^'^ - 2^sme cos{(p - wti)y^^ . (4.5) 

Similarly, since the time interval t — t2 that needs 62 to go from its retarded position 
C2 to the actual position A2, is the same that takes the light for travel from C2 to 

P, we have 

t-t2 = ia/c)C^[l + f + 2^smecosi(p-ujt2)y^\ (4.6) 
where ^ denotes the parameter 

^ = - < 1. (4.7) 
r 

Equations (4.5) and (4.6) are complicated functional equations that determine in an 
unique way the retarded times ti and t2 respectively, as functions of the parameters 

t, r, 6 and (/?. The retarded times ti and t2 are the same only for detection points 
over the z axis. Instead of working with the retarded times ti and t2, it is convenient 
to introduce the following variables. 

X = (f — uti, (4-8) 

y = (p - ut2. (4.9) 

In the integral (4.1) the parameters t, r and 9 are fixed; then Eq. (4.5) determines 
ti as a function of the angle Lp, and therefore the variable x defined in Eq. (4.8) has 
a unique value for each p in the interval < < 27r. This property allows us to 
carry out the integral (4.1) as an integral over the variable x. In order to see this 
clearly, let us first note that the correspondence between the variables x and y in 
Eqs. (4.8), and (4.9) is one to one. In fact, from Eqs. (4.5) and (4.6) we get the 
following relation between x and y 

y - X = I (1 + + 2^ sin 61 cos j/) - (l + - 2^ sin 9 cos | (4.10) 

where f3 = au/c. Taking the derivative with respect to x in Eq. (4.10), we obtain 

dy 1 + pj~^/3sin6'sinx ("4 11) 

dx 1 + p^^l3sm.9smy' 

where 

= (1 + ^2 _ 2^sin6/cosa;)^/^ (4.12) 
9 



and 



But 



and 



therefore 



P2 = (l+^^ + 2esin^cosj/)^/2 (4 13) 



1 + Psm9smx = Ki > 0, (4.14) 



1 + /92 ^/3sin6'siny = K2 > 0, (4.15) 



^ = =i>0, (4.16) 

ax K2 



which proves that the correspondence between x and y is one to one. This property 
holds for any time t, radius r and angle 9. 

Let us consider now the integral (4.1); where the time t, the radius r and the 
angle 6 remain fixed. From Eq. (4.5) we obtain 

-u;^ = (p^'(3sm9smx) ^; (4.17) 
dip ^ ^ dip 

if we combine this equation with Eq. (4.8), we get 

dx 1 / . X 

— = — > 0, 4.18 
dip K\ 

where Ki is explicitly given in Eq. (4.14). Now, the equation dip/dx = ki > 
tells us that ip is an strictly monotonous increasing function of x; so we can put the 
integral (4.1) in the following form 

27r+a 

j (El X B2 + E2 X Bi) • TKidx, (4.19) 

a 

where the parameter a is given by 

a = -uti{ip = 0) = -i^ti{ip = 2Tr), (4.20) 

which in general depends in a complicated way on the time t, the radius r and the 
angle 9 of the band. When the integrand of (4.19) is explicitly evaluated by using 
the electric field (4.2) and the magnetic induction (4.3), with the corresponding 
expression for E2 and B2, it can be shown that the variables x and y appear only 
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as sinx, cosx, siny and cosy. The integrand is, of course, a function of the 
variable x only since, as shown above, y is uniquely determined by the value of x in 
Eq. (4.10). Moreover, as it can be easily seen, the correspondence between x and y 
is such that if y is the value associated with x, then y + 27r is the value associated 
with X + 2tt. Thus we conclude that the integrand of (4.19) is a periodic function of 
x, with a period of 27r. This property implies at once that the integral (4.19) does 
not depend on the value of the parameter a; so we can put a = in it, obtaining 
the following representation for the energy- flux (4.1). 

27r 

j (El X B2 + E2 X Bi) • TKidx. (4.21) 


In particular then, the energy flux across the band between 9 and 6 + dO over the 
sphere of radius r docs not depend on the value of the time at which it is evaluated. 
The integral (4.21) depends, however, in a very complicated way on the radius r 
and the angle 9. 

The time-independence of the integral (4.21) is true for any band over the surface 
of the sphere of radius r; therefore the energy flux across the whole surface of the 
sphere is also independent of time. This property was inferred on symmetry grounds 
in section II. There we also proved that the interference of radiation given by 

7r 27r 

(c/47r)r2 j sm9d9 j (Ei x B2 + E2 x Bi) • rKidx (4.22) 


can be evaluated for an arbitrary radius r with r > a, because it is independent of 
r. The integrand of (4.22) contains a great number of terms for any finite value of r, 
and in addition Eq. (4.10) that links the variables x and y is very complicated for an 
arbitrary r. Strong simplifications of the integrand of (4.22) and of the functional 
relation (4.10) are obtained when considering the limit when r goes to infinity. We 
emphasize the fact that, due to the independence of the interference rate on the 
radius r, this limit does not present any complication, being perfectly well defined. 
In this limit Eq. (4.10) is reduced to 

y — x = (3sm9{cosy + cosx), (4.23) 
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and the interference rate of radiation (4.22) becomes 



(4.24) 

with 7(/3) given by 

7r/2 2n 

3 f „ ,„ / (cos^ ^cosxcos-u + sinxsiny — sin^ 6)(ia; 
47r J i (1 - /3sin6'sinx)2(l + /3sin6'siny)3 ^ ' 

If instead of changing the variable ip in the integral (4.1) by the x of Eq. (4.8), 
we perform the integral (4.1) by means of the variable y defined in Eq. (4.9), we 
would obtain in place of (4.25) the following expression for I {13). 

3 f . f {cos^d COS X COS y + sin xsiny- P^sin^ e)dy 

'^^^ = ^ I ""^^^ / (l-/3singsinx)3(l+/3singsiny)^ " ^^'^'^ 

This formula looks different from (4.25), but it can be easily shown, with the help of 
Eq. (4.23), that both formulae are the same. We also point out that the functional 
Eq. (4.23) has the solutions 

x = y = 7r/2 (4.27) 

and 

x = y = 37r/2 (4.28) 

independently of the value of /? and 9. 

Eq. (4.25), with y(x) defined in Eq. (4.23), is an exact formula for the inter- 
ference radiation, and it is the main result of this paper. Unfortunately, due to the 
complicated relation between the variables x and y defined by Eq. (4.23), the inte- 
gral (4.25) cannot be evaluated in a closed analytical way. This situation contrasts 
with what happens when we consider the case of one charge in circular motion using 
the present treatment, where the corresponding integral can be explicitly done as it 
is shown in appendix B. 

The integral (4.25) can be easily studied by means of numerical techniques. In 
Fig. 3 we show the function y{x) in the orbit plane, for three different values of the 
parameter /3; namely for /? = 0, /3 = 0.5 and /3 = 0.9999. In Fig. 4 we present the 
results of the numerical treatment of the function /(/?) of Eq. (4.25) in the interval 
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< /? < 0.99 [14|, where instead of the variable /? we have used the variable (j) 
defined by 

/3cos(/> = (4.29) 



since it is more convenient. From the last equation we obtain 

d(3 1+I3sui4> 



(4.30) 



dcp cos (p 

The condition j3 < 1 implies cj) < 0.739, so that dp/dcp > 0. This shows that the 
correspondence between (3 and (j) is one to one. 

In figure 4 we have also drawn a dotted curve that represents the power series 
ofEq. (3.5) 



FIGURE 3 



FIGURE 4 



V The interference rate in the ultrarelativis- 
tic case 

When f3 is very close to one, the integrand of (4.25) is significative only for values 
of the variables x and y around 

XI = yi = 7r/2 (5.1) 



X2 = y2 = 3^2 (5.2) 
In the approximate evaluation of (4.25) we will use a procedure similar to that of 



reference [15[ for the one charge case; but now the approximations are more crude 
because the integrand of (4.25) around (5.1) and (5.2) is not as sharply defined as 
in the one electron case. Nevertheless, this somewhat heuristic procedure allows us 
to obtain a simple analytical formula for the leading part of the interference rate, 
which accuracy improves according as (3 becomes close to one. 

Since the radiation is mainly concentrated in the orbit plane, it is convenient to 
introduce the angle 

X = 7t/2-9 (5.3) 
13 



Then, we are going to consider an expansion of the integrand of (4.25) with dx = 
X — xi and x of the order of 7"^. Prom (4.23) it fohows then that dy = y — yi is of 
the order of 7"^, that is, y practically does not change when x is around 7r/2, and 
we can put y = 7r/2. In this way we can approximate the integral of x around xi in 
(4.25) by 

Y'^ [(l/2)(x-xi)2-a2]dx 

J 2[(x-xi)2 + a2]2 l^-^^ 

xi—Sx 

where 

= 7-'(l + bx?). (5.5) 

Extending the limits of integration between —00 and +00 in (5.4), we obtain for it 
the value it /8a. Now, since the radiation is mainly concentrated in the orbit plane, 
we can approximate sin^ by 1 in the outermost integration of eq (4.25), and if the 
variable x in changed by z = 7%, we get the following contribution around xi = tt/2, 
for the interference rate of radiation. 

77r/2 


The contribution l2iP) around X2 = 37r/2 of the integral (4.25) must be, on symme- 
try grounds, equal to (5.6). In this way, we get the following approximated formula 
for the interference rate of radiation in the ultrarelativistic case. 

- ^^-^w 

where 

/(7) = ^ln(77r) (5.8) 

In table 5.1 we represent the interference rate given by the approximated formula 
(5.8), and the value of this quantity evaluated numerically from the exact formula 
(4.25). As expected, the accuracy of (5.8) improves when 7 increases. 
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Table 5.1 


7 




/(7) 


Accuracy 


100 


1,235 


1.0781 


-12,7 % 


500 


1,536 


1,3799 


-10,16 % 


1000 


1,666 


1,5098 


-9,37 % 


5000 


1,968 


1,8116 


-7,95 % 


10000 


2.098 


1,9416 


-7.45 % 


30000 


2.301 


2,1170 


-6.79 % 



The behavior (5.8) for the interference rate of radiation in the ultrarelativistic 
case contrasts strongly with the behavior of the Larmor terms of each charge, which 
behaves as 7^ when 7 goes to infinity. In particular, in the ultrarelativistic case, 
the interference rate of radiation is completely negligible in comparison with the 
Larmor terms of each charge. 
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APPENDIX A 

The far field produced by two charges in arbitrary motion can be expressed in the 
following form: 

oo 

E = ^ E„ (^.1) 

n=l 

where E„ is the last term in Eq. (3.1). Introducing the notation 



r = ri = — r2 = a cos uti + a sin utj (A. 2) 

V = Vi = —V2 = —aoj sin UJU + au cos ojtj (A. 3) 

a = r ■ r = ar sin 9 cos{ip — oot) (A.4) 

rj = r -v = aujr sin 6 sin{ip — ut) (A-.5) 

r = sin cos tpi + sin 9 sin ip hatj + cos 9k (A-.6) 



we obtain the following expressions for the fields: 



= ^^^^{r-af} (A.7) 

E2 = ^^^1±|^ {r/r + av - 2a?7f } (A.8) 

Es = {(6??^ - 7uj^a^)r + Uariv + {7u^a^ - 20ar?^)r} (A.9) 

E, = ^^^^ + J^^' {(37?=^ - 12u;^a%)r + (12r?^a - 5a;^a^)v 

+ {ULu'^a^r] - 15a7?^)f } (A.IO) 

E, = i^i-Lf^|(i20r?4_i080a;Va' + 241a;V)r+ 
(720?7^a - 964oJ'^a^r])^r+ 

{-SAOri'^a + 2{)A^'^rj^a^ - 241a;^a^)f } (A.ll) 

Eq = l£l±^4^((720r/5- 12000a;2aV + 8560u;4Q4r?)r+ 

120rc' 

(eOOOar?'^ - 17120a;2a^r/2 + 2256u;^a^)v + 

(-6720ar?^ + 29120uj'^ a^r]^ - 10816a;^a^r/)f } (A.12) 
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El = ( (5040a6 - 138600a;2 alpha^rf + 209790a;'' a^'r/^ - 19279a;^a^)r+ 

720rc* 

(55440ar/^ - 2797200;^ a^r/^ + limi Aw'^a^r])^f + 

(-60480ary^ + 418320a;2a^r?^ - 325464 omega^a^if + 19279a;^a^)f } (A. 13) 

= M%±^^|(_ioi50a;6aS + 35280a;4aV- 13230a;2a2^5^3157?^)r+ 
315rc^ 



(-1957a;^a^ + 30450a;^a^7?2 - 35280a;2a^r?^ + 4410ar/'')v + 



(12107w^a^r/ - 65730w^a^r/^ + 48510a;^a^r/^ - 4725a7?'^)f } (A.14) 

Eg = '-^^ ~ ^^^^^ I (27715210;^ - 55597920^^0^2 + 92786400a;^a^r/^- 
40320rc^^ 



22014720a;2a%^ + 362880?7^)r + 

(-22172168a;^a^r? + 111195840a;'' a^T?^ _ 742291200;^ a=^r?^ + 6289920ar?^)v 
(-2771521a;^a^ + 77770088a;^a^r/2 - 203982240a;^aS^ + 96243840a;2a3r?^ 

-66528000;??'^ )f} (A.15) 

For circular motion the Pointying vector becomes 

S = -^EH M.15) 
47r 

and it can be evalued using the expressions of this appendix. If we now perform an 
integration over the surface of a sphere centered in the origin and that encloses the 
orbit, we obtain eq. (3.3). 
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APPENDIX B 



In the case of one electron in circular orbit with constant velocity, the energy flux 
across the spherical surface or radius r and center at the orbit center, namely 

(crV47r) J sm9de J (Ei x Bi) • rdip, {B.l) 



is, like Eq. (2.4), independent of the time t and of the radius r of the sphere. On 
introducing in (B.l) the variable x of Eq. (4.8), and on taking the limit when r 
goes to infinity, it becomes 

/ 2/. 2Na4 ? • /.,^Tf 1 sin^ cos^ X 1 , 

(ce^/47ra^)/3* / smOdO / { ^ . ^ . - -ji — ^ ■ n ■ t dx 

^ ' ' J I (1 - /3sin6'sina;)3 (1 - /3sin6lsinx)5 J 

These integrals, unlike the integrals of Eq. (4.25), can be easily evaluated in a closed 
analytical way. The result is, of course, Eq. (3.4). 
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Figure and Table Captions 



Table 5.1. Comparison between the numerical calculation of the exact expression 
(4.25), here denoted by 1(7), and the asymptotic formula (5.8), denoted 
by 1(7). We also show the porcentual error with respect to the numeri- 
cal value. As expected, the approximation (5.8) improves as 7 becomes 
larger. 

Figure 1. Two charges ei and 62 in circular motion at constant angular velocity 
ijj. The orbit has radius a, and the coordinate axes have been chosen 
so that the orbit is centered at the origin and contained in the X — Y 
plane. We also show two spherical surfaces Si and S2 centered at the 
orbit's center, O being the domain bounded by Si and S2. 

Figure 2. At the observation, time t the charges e\ and 62 are ubicated at opposite 
ends of diameter A1A2. In the same way, at the retarded times ti and ^2 
of charges ei and 62 they are ubicated at the ends of the dotted diame- 
ter B1B2 and the dashed diameter C1C2 respectively. We also show the 
retarded distances Ri and R2 from the retarded positions of the charges 
to the observation point P, and the radius r of the spherical surface to 
which P belongs. 

Figure 3. The curves represent the function y{x) defined by the retardation con- 
dition, eq. (4.23). The plot is made for three different values of the 
parameter namely P = 0, /3 = 0.5 and /9 = 0.9999. Note that in all 
cases X = y = 7r/2 and x = y = 3tt/2 are solutions of eq. (4.23). 

Figure 4. The solid lines show the numerical evaluation of function /(/?) defined 
by eq. (4.25). The plot is made in terms of the variable (f) defined by 
13 cos (j) = (j), for the range < /3 < 0.99. We also show a dotted line 
that represents the power series expansion for the interference rate, eq. 
(3.5), in the range < /3 < 0.4. 
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